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Using a sample of 383 million BB events collected by the BABAR experiment, we measure sums of 
seven exclusive final states B -^ Xd(s)J, where Xd{Xs) is a non-strange (strange) charmless hadronic 
system in the mass range 0.6 — 1.8GeV/c^. After correcting for unmeasured decay modes in this 
mass range, we obtain a branching fraction for & — > 1^7 of (7.2 ± 2.7{stat.) ± 2.3{syst.)) x 10^®. 
Taking the ratio of Xa to X^ we find r(& -^ d'y)/T{b ^ S7) = 0.033 ± 0.013(stat.) ± 0.009(syst.), 
from which we determine \Vtd/Vu\ = 0.177 ± 0.043. 

PACS numbers: 13.20.Hc 



The decays b ^ d'y and b —^ sj are flavor-changing 
neutral current processes forbidden at tree level in the 
standard model (SM). The leading-order processes are 
one-loop electroweak penguin diagrams in which the top 
quark is the dominant virtual particle. In the SM the in- 
clusive rate for b ~^ dj is suppressed compared to 6 ^ sj 
by a factor of iVtrf/Visp, where Vtd and Vts are Cabbibo- 
Kobayashi-Maskawa matrix elements. Measurements of 
|Vtrf/Vti| using the exclusive modes B —f {p,Lu)y and 
B -^ K*j [1] have theoretical uncertainties of 7% from 
weak annihilation and hadronic form factors [2] . A mea- 
surement of the inclusive decay 5 — > (i7 relative to 5 — > 57 
could determine |Vtd/Vts| with reduced theoretical uncer- 
tainties compared to the exclusive modes [3[ . In theories 
beyond the SM 4] , new virtual particles may appear dif- 
ferently in the penguin loop diagrams for b —> dj and 
b ^ sj and in the box diagrams responsible for B^ 
and Bg mixing [5'], leading to measurable differences in 
|Vtd/Vts| extracted from these two methods. 

We present measurements of the rare decays B -^ Xdj 
using seven exclusive final states (see Table U) in the 
hadronic mass range 0.6 < M{Xd) < l.OGeV/c^ (which 
contains the p and lu resonances), and in the previously 
unmeasured region 1.0 < M{Xd) < 1.8GeV/c^. We 
combine our results in the two mass regions and make 
corrections for decay modes that are not reconstructed 
to obtain an inclusive branching fraction for b —f d'y 
in the mass range 0.6-1.8 GeV/c^. We perform a par- 
allel analysis of _B — > Xs7 using the equivalent seven 
modes (Table IJ), and determine the ratio of inclusive 
rates T{b ^r d'y)/T{b -^ sj) in the hadronic mass range 
0.6< M(Xrf) < 1.8GeV/c2. 

These measurements use a sample of 383 x 10^ BB 



TABLE L The reconstructed decay modes. Charge conjugate 
states are implied throughout this paper. 



B^Xdj 




B^Xs-f 


B" ^ 7r+7r-7 
B+ ^ 7r+^°7 
B+ ^7r+7r"7r+7 

B° ^ TT + ^-TT^ 
B+ ^ TT + ^-TT + V 

B+ -^ 7r+777 


'7 
'7 


B° ^ K+TT-'y 
B+ ^ K+TT^-y 
B+ ^ K+TT-Tv+'y 

B" ^ K+7T-TV°J 

B° -* K+n-7v+n-y 

B+ -^ K+TT-TV+TT°J 

B+ ^ K+rj-y 
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pairs collected at the T(4S') resonance with the BaBAR 
detector Q at the PEP-II B factory. The high-energy 
photon is reconstructed from an isolated energy cluster 
in the CsI(Tl) calorimeter, which has a shape consistent 
with a single photon, and an energy 1.15 < E* < 3.5 GeV 
in the center-of-mass (CM) frame. We remove photons 
that can form a tt^ (77) candidate with another photon 
of energy greater than 30(250) MeV, if the two-photon 
invariant mass is in the range 105 < m^^ < 155 MeV/c^ 
(500 < m^^ < 590MeV/c2). 

Charged pion and kaon candidates are measured in a 
1.5 T magnetic field as tracks in a 5-layer silicon vertex 
detector and a 40-layer drift chamber, with a minimum 
momentum in the laboratory frame of 300 MeV/c. To dif- 
ferentiate pions from kaons we combine information from 
the detector of internally reflected Cherenkov light with 
the energy loss measured in the tracking system. At a 
typical pion energy of 1 GeV, the pion selection efficiency 
is 85% and the kaon mis-identification rate is 3%. Kaons 
are selected by inverting the pion selection criteria. We 
reconstruct 7r°(r/) candidates with momenta greater than 
300 MeV/c from pairs of photons of minimum energy 20 
MeV with an invariant mass 107 < m^-y < 145MeV/c^ 
(470 < m^j < 620MeV/c^). The selected charged tracks, 
T^^ijf) candidates, and high-energy photons are combined 
to form B meson candidates consistent with one of the 
seven B -^ Xgj 01 B ^ Xd'y decay modes. For B -^ Xg^ 
decays one charged kaon is required, with all other tracks 
required to be pions. For B — > X^j decays, all tracks are 
required to be identified as pions. The charged particles 
are combined to form a common vertex with a vertex fit 
probability greater than 2%. 

Most of the backgrounds in this analysis arise from 
continuum e+e^ —> qq events, q — (u, d, s, c), in which a 
high-energy photon comes from either initial state radia- 
tion or the decay of a TT^irj) meson. We require R2 < 0.9 
and I cos 6*7 1 < 0.8, where R2 is tjbe ratio of the second 
to zeroth Fox- Wolfram moments [7] , and 9t is the angle 
between the photon and the thrust axis of the rest of the 
event (ROE) in the CM frame. The ROE includes all the 
charged tracks and neutral energy in the calorimeter not 
used to reconstruct the B candidate. 

The quantity cos 9t and twelve other variables that dis- 
tinguish between signal and continuum events are com- 
bined in a neural network (NN). These include the ra- 
tio i?2, which is i?2 is calculated in the frame recoiling 
against the photon momentum, the B meson production 
angle 9g in the CM frame with respect to the beam axis, 
and five Legendre moments of the ROE with respect to 
both the thrust axis of the ROE, and the direction of 
the high-energy photon. Differences in lepton and kaon 
production between background and B decays are ex- 
ploited by including five flavor-tagging variables applied 
to the ROE 8|. We optimize the NN configuration for 
maximal discrimination between signal and continuum 
background, which gives 50% signal efficiency and 0.5% 
misidentification of continuum background. 



We use the kinematic variables AE = Et 



-^bcamand 



?*2 



where Eg and p^ are the CM 



energy and momentum of the B candidate, and £'^(,^111 
is the CM energy of one beam. Signal events are ex- 
pected to have a /S.E distribution centered at zero with 
a resolution of about 30 MeV, and an tties distribution 
centered at the mass of the B meson with a resolution of 
about 3 MeV/c^. We consider candidates in the ranges 
-0.3 GeV < A£ < 0.2 GeV and wes > 5.22GeV/c2 to 
incorporate sidebands that allow the combinatorial back- 
ground yields to be extracted from a fit to the data. On 
average there are 1.75 candidates per event, and in events 
with multiple candidates we select the one with the best 
TT^{r]) mass, or, where there is no iT^{rj) we select the 
candidate with the best vertex fit probability. 

The signal yields in the data are determined from two- 
dimensional unbinned maximum likelihood fits to the A_E 
and rriES distributions of the sums of all seven final states 
listed in Table HI We consider the following contribu- 
tions: signal, combinatorial backgrounds from continuum 
processes, B — > Xn^ /r] decays, backgrounds from other 
B decays, and cross-feed from mis-reconstructed signal 
B — > X7 decays. The fits to the B -^ X^^ samples con- 
tain a component from misidentified B — > Xs^ decays, 
but we neglect the small B -^ Xdj background in the 
B — > Xs^ samples. The B background yields are deter- 
mined from a Monte Carlo (MC) simulation, whereas the 
continuum background is allowed to float in the fit. 

Each background contribution is modeled by a prob- 
ability density function (PDF) that is determined from 
MC. The signal PDFs are the product of one-dimensional 
rriES and AE distributions determined from fits to the 
B -^ K*^ data. For the signal cross-feed component, 
and the B -^ Xgj background in the B -^ Xd'f fit, 
we use two-dimensional histogram PDFs to account for 
correlations. The contributions from B —i- Xtt^/i] are 
modeled by Gaussian peaks in both AE and toeS: where 
AE is displaced by —80 MeV due to the missing photon. 
The B — > Xs7 background in the B -^ Xd'y sample also 
peaks with AE displaced by —50 MeV due to the kaon 
misidentification. Continuum and other non-peaking 
backgrounds are described by an ARGUS shape [9] in 
Toes and a second-order polynomial in AE. 

We perform fits separately for B — » X^j and B -^ Xgj 
and in the two hadronic mass ranges. The signal and con- 
tinuum yields and the ARGUS and polynomial contin- 
uum shape parameters are allowed to vary. We scale the 
cross-feed contribution proportionally to the fitted signal 
yield, re-fit, and iterate until the fit converges. The fits 
for B —> Xsj and B —>■ Xd'y are shown in the low- and 
high-mass regions in Figs. [T] and [21 respectively. 

The signal yields, average efficiencies and partial 
branching fractions for the sums of the seven decay modes 
are given in Table [TTl The reconstruction efficiency de- 
pends on the distribution of the signal yield among the 
final states. For Xg we measure the distribution of the fi- 
nal states in the data, but for Xd there is no statistically 
useful information, so we model the distribution using 
the the phase space fragmentation model implemented 




FIG. 1: Projections of the fits to data in the hadronic mass 
range 0.6-1.0 GeV/ c^. Projection of AE with 5.275 < rriES < 
5.286 GeV/c^ for (a) B -^ Xsj and (c) B -^ XdJ, and mES 
with -0.1 < AE < 0.05 GeV for (b) B -^ Xs-y and (d) B ~> 
Xd'y- Data points are compared with the sum of all the fit 
contributions (solid line) including the signal (dashed line). 
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FIG. 2: Projections of the fits to data in the hadronic mass 
range 1.0-1.8 GeV/ c^. Projection of AE with 5.275 < tties < 
5.286 GeV/c^ for (a) B -^ Xs-y and (c) B ~+ Xd'y, and jties 
with -0.1 < AE < 0.05 GeV for (b) B -^ Xsj and (d) B -^ 
Xd'y- Data points are compared with the sum of all the fit 
contributions (solid line) including the signal (dashed line). 



in JETSET ^. 

The branching fractions in Table IIIII are obtained af- 
ter correcting for missing final states. The low mass 
B — > Xs"f measurement is found to be consistent with 
previous measurements of the rate for B -^ K*j [lll |. 
after accounting for the 50% of decays to neutral kaons. 
For the low mass B -^ Xd'y region, non-reconstructed p 
and to decays are small and we find a branching fraction 
of (1.2 ± 0.5) X 10~^, consistent with previous measure- 
ments of B{B -^ {p,uj)^) [l|. In the high mass region, 
we correct for missing final states with > 5 stable parti- 
cles, or with multiple tt'^s, using the fragmentation model 



TABLE II: Signal yields (Ns), average efficiencies (e) and 
partial branching fractions (B) for the measured decay modes. 
The first error is statistical, the second systematic. 



M(X)[GeV/c 



Ns 



B(xl0-'') 



0.6 < M{Xs) 
0.6 < M{Xd) 
1.0 < M{X,) 
1.0 < MjXd) 



< 1.0 

< 1.0 

< 1.8 

< 1.8 



1543 ± 46 

66 ±26 

2279 ± 75 

107 ± 47 



8.5% 
7.0% 
6.1% 
5.2% 



23.7 ±0.7 ±1.7 
1.2 ±0.5 ±0.1 

48.7 ±1.6 ±4.1 
2.7 ±1.2 ±0.4 



TABLE III: Branching fractions B(xlO"'') in the two 
hadronic mass regions M(X)[GeV/c'^], after correcting for 
missing final states, and the ratios of B{b —* d-y) to B{b —* 57) . 
The first errors are statistical, and the second are systematic, 
including the fragmentation of the hadronic system. 



M{X) 


B{h ^ d'y) 


B{b -^ s-f) 


B{b -^ d-f)/B{b -^ s-f) 


0.6- 1.0 
1.0- 1.8 
0.6-1.8 


1.2 ±0.5 ±0.1 
6.0 ±2.6 ±2.3 
7.2 ±2.7 ±2.3 


47 ± 1 ± 3 
168 ± 14 ± 33 
215 ± 14 ± 33 


0.026 ±0.011 ±0.002 
0.036 ±0.015 ±0.009 
0.033 ±0.013 ±0.009 



described above. 

The sources of systematic uncertainties in the measure- 
ment of the branching fractions are listed in Table IIVI 
These include uncertainties on track reconstruction effi- 
ciency, 7 and 7r°/?7 reconstruction, the tt^/t; veto, the NN 
selection, and the number of BB pairs. The 2% uncer- 
tainty on correct kaon/pion particle identification, and 
the 20% uncertainty on kaon misidentification, which is 
a systematic on the fixed b ~> sj background in the 
B -^ Xij'j fits, do not cancel in the ratio. The system- 
atic errors associated with the variation of the fit PDFs 
also do not cancel because of the very different signal 



TABLE IV: Systematic errors on the measured partial and 
total branching fractions B. The final column shows system- 
atic errors that do not cancel in the ratio of rates r(6 — > 

d7)/r(&->s7). 



Systematic 


M{Xs) 


M(Xd) 


Xd/Xs 


Error Source 


0.6-1.0 


1.0-1.8 


0.6-1.0 


1.0-1.8 


Ratio 


Tracking 


1.7% 


1.7% 


1.7% 


1.7% 




High-energy photon 


2.5% 


2.5% 


2.5% 


2.5% 




7r°/?7 reconstruction 


1.7% 


1.7% 


1.7% 


1.7% 




7r°/ry veto 


1.0% 


1.0% 


1.0% 


1.0% 




K/tv identification 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


Neural network 


5.0% 


5.0% 


5.0% 


5.0% 




BB pair counting 


1.1% 


1.1% 


1.1% 


1.1% 




Fit PDFs 


2.4% 


3.6% 


7.0% 


8.3% 


8.7% 


Backgrounds 


0.3% 


0.4% 


2.4% 


6.1% 


5.4% 


Fit bias 


0.4% 


1.7% 


0.4% 


3.3% 


3.0% 


Fragmentation 




3.6% 




7.7% 


8.5% 


Partial B 


7.0% 


11.4% 


10.0% 


14.8% 


13.8% 


Missing > 5 body 




5.6% 




25.8% 


21.0% 


Other missing states 




17.0% 




23.8% 


7.1% 


Spectrum Model 




1.8% 




1.6% 




Total B 


7.0% 


21.2% 


10.0% 


38.1% 


26.1% 



to background ratios in the two samples. We vary the 
signal PDF parameters within the range allowed by the 
fit to the B -^ K*^ data. The normalization of the sig- 
nal cross-feed is varied by ±30%, and the contribution 
of B ^ Xtt^/tj by ±100%, in accordance with MC stud- 
ies. The remaining peaking B backgrounds, including the 
B —> Xs^ contribution in the B — > ^^7 sample, are var- 
ied by ±20%. We use simulated signal and background 
event samples to assign a systematic uncertainty due to 
the potential for bias in the fit method. 

There is an additional systematic error on the efficiency 
due to the uncertainties in the measured fragmentation 
of the Xs hadronic system into the seven B — > Xg^i fi- 
nal states. The equivalent error for B -^ Xj^j is obtained 
from the difference between our fragmentation model ap- 
plied to B ^ Xdj and the fragmentation observed in 
B -^ Xsj data. We assume that these errors are inde- 
pendent and so do not cancel in the ratio of branching 
fractions. 

Table IIVI also shows the systematic errors associated 
with correcting the partial branching fractions for the 
missing final states. There is no information from the 
data on the missing fraction of high multiplicity final 
states with > 5 stable hadrons, or on the missing fraction 
of other final states with > 1 7r° or 77 mesons. We vary 
these fractions by ±50% of their values from the default 
phase space fragmentation. We motivate our choice of a 
±50% variation using signal models, for which we mix a 
combination of resonances as 50% fractions of S ^ Xgj 
and B — > Xd'y in the mass range 1.0 — 1.8 GeV/c^. These 
give missing fractions close to the lower limits from the 
±50% variations. The missing fraction errors partially 
cancel in the ratio when the ±50% variations are made 
in the same direction for b —f d'y and b ~^ sj. 

We take the spectral shape of the high-energy 
photon from I12l| with the kinetic parameters 
(TOfc,^2) ^ ' (4.65 GeV/c^, -0.52 GeV^) extracted 
from fits to 6 — > 57 and b — > c£iy data [13|. We 
vary these shape parameters in a correlated way 
between (mb,/!^) = (4.60 GeV/c^, -0.60 GeV^) and 
{mb,fil) = (4.70 GeV/c2, -0.45 GeV^). There are 
systematic errors on the branching fractions from 
these variations, but they are small and cancel in 
the ratio. The fraction of the spectrum in the mass 
range 0.6-1.8 GeV/c^ 
(50 ± 4)% for b -^ 
ratio of branching fractions to Mx > 1.8GeV/c^, so 
these errors, which mostly cancel in the ratio, are not 
included in Table IIVI If we make this correction, we 
obtain Bib -^ d'y) = (1.4 ± 0.5 ± 0.4 ± 0.1) x 10"^ and 
B{b -^ 57) = (4.3 ±0.3 ±0.7 ±0.2) x lO"*, where the first 
error is statistical, the second systematic and the third 
accounts for the uncertainty in extrapolating to the full 
mass range. The result for B — > Xgj is consistent with 
the measured inclusive b —> sj branching fraction of 
(3.55 ±0.24) X 10-^ ;il]. 

We convert the ratio of branching fractions from the 
full mass range 0.6-1.8 GeV/c^, F(6 -> d'y)/T{b -> 57) = 



is (51 ± 4)% for b ^ d'y and 
sj. We do not extrapolate the 



0.033 ± 0.013 ± 0.009, into a value for \Vtd/Vts\ us- 
ing Table 1 and Equation (26) of [3]. The result is 
\Vtd/Vts\ = 0.177 ± 0.043 ± 0.001, where the first error 
is experimental, including systematic errors, and the sec- 
ond error is theoretical. The theoretical error includes 
uncertainties on the CKM parameters p and 77, and on 
1/to^ and I/ttt,^ corrections, but does not include an un- 
certainty for the restriction of the measurement of the 
ratio to hadronic masses below 1.8 GeV/c^. 

As a check, we use the low mass region to determine 
I^td/Vtsl using predictions for exclusive B — > {p,u})^ and 
B -^ K*j from 0. We find \Vtd/Vts\ = 0.214 ± 0.046 ± 
0.028 where the first error is experimental and the second 
is theoretical. This is in good agreement with previously 
published results [l|. 

In summary we have made the first measurement of 
B -^ Xdj decays in the hadronic mass range up to 



1.8GeV/c^, and have extracted |14d/^ts| from an inclu- 
sive model with small theoretical uncertainties. These 
results are consistent with the measurements of |14d/Vts| 
from the exclusive decays B -^ {p,uj)^ [1], and with 
Bs/Bd oscillations Q- 
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